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Dissipation and decoherence, and the evolution from pure to mixed states in quantum physics are 
handled through master equations for the density matrix. Master equations such as the Lindblad 
equation preserve the trace of this matrix. Viewing them as first-order time-dependent operator 
equations for the elements of the density matrix, a unitary integration procedure can be adapted to 
solve for these elements. A simple model for decoherence preserves the hermiticity of the density 
matrix. A single, classical Riccati equation is the only one requiring numerical handling to obtain a 
full solution of the quantum evolution. The procedure is general, valid for any number of levels, but 
is illustrated here for a three-level system with two driving fields. For various choices of the initial 
state, we study the evolution of the system as a function of the amplitudes, relative frequencies and 
phases of the driven fields, and of the strength of the decoherence. The monotonic growth of the 
entropy is followed as the system evolves from a pure to a mixed state. An example is provided by 
the 71 = 3 states of the hydrogen atom in a time-dependent electric field, such degenerate manifolds 
affording an analytical solution. 

PACS numbers; 03.65.Yz, 05.30.-d, 42.50.Lc, 32.80.Qk 



UNITARY INTEGRATION PROCEDURE FOR 
MASTER EQUATIONS 

Master equations, such as the Lindblad equation 0, 
can describe dissipation and decoherence in quantum sys- 
tems. In recent work [3^, one of us adapted a "unitary 
integration" procedure [3, Q for solving such equations 
while preserving desirable properties such as the her- 
miticity of the density matrix even in the presence of 
dissipation and decoherence. This permits keeping track 
of quantities such as the entropy while the system evolves 
from a possibly initial pure state to a final mixed one. 
The two-state illustration given in that initial work is 
extended now to a three-level system through suitable 
combinations of density matrix elements to preserve the 
hermiticity of the operators involved. 

Consider the master equation for the density matrix 
p called the Liouville-von Neumann-Lindblad equation 

ip = [H, p] + [[Lkp, lI] + [ifc, pill) 

k 

= [H, p] - {^i^f^P + P^l^^ - ^LkpLi) ,(1) 

k 

where an over-dot denotes differentiation with respect to 
time and h has been set equal to unity. H is a Hermi- 
tian Hamiltonian while the Lj, are operators in the sys- 
tem through which dissipation and decoherence are intro- 
duced. Even though this can result in non-unitary evo- 
lution, the form of the equation preserves Tr(/5) and pos- 
itivity of probabilities. A more mathematical discussion 
of such "super-operators" and "dynamical semigroups" 
is given in 



A commonly used form of H is 

Hit) = e{t)A, + 2J{t)A,, (2) 

with 
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The couplings indicated in Eq. Q between states 1 and 
2 and between 2 and 3 of a three-state system are re- 
ferred to as A and V depending on the relative energy 
positions of the three states, whether 2 lies above or be- 
low, respectively, relative to levels 1 and 3. The three 
operators in Eq. Q close under commutation according 
to the standard relations satisfied by angular momentum 
algebra: [^a,, A^] = iA^, and cyclic. Hioc and Eberly 0| 
considered such a Hamiltonian for the Liouville version of 
Eq. 1^, that is, without the dissipative term, along with 
solutions for certain forms of e and J. Population trap- 
ping and dispersion was also considered in with a sim- 
ilar Hamiltonian, and |^ generalized to n-level systems. 
Our work presented here may be regarded as extending 
such studies to include also dissipation and decoherence. 

In general, with each of the three states having distinct 
energies Ei, E2, E^, and the driving fields having finite 
detunings from resonance, entries along the diagonal of 
H in Eq. Q complete the Hamiltonian for such systems. 
Full treatment according to our formalism below then re- 
quires all the elements of the SU(3) algebra, namely five 
more linearly independent 3x3 matrices to supplement 



2 



those in Eq. (jj)). We expect to return to this later but, 
in this paper, we restrict ourselves to the degenerate case 
of equal eigenvalues in which case the above three matri- 
ces suffice and the calculations reduce to solving a single 
equation just as in the two-state system considered in 
01 • Applications include three identical coupled pendula 
with nearest neighbor time-dependent couplings, driven 
systems on resonance, and the degenerate states of the 
n — S manifold of hydrogen driven by time-dependent 
electric fields. 

Dissipation and decoherence are introduced through 
the Lk matrices in Eq. Here again, as shown in 2], a 
choice of all eight linearly independent matrices affords a 
simplification because of a sum rule that inserts the de- 
coherence as a unit operator in such an eight-dimensional 
space. In this procedure, Eq. ^ is recast into a set of 
eight equations for the elements of the density matrix 
(recall that the trace remains invariant). An appropriate 
linear combination of the elements such that the opera- 
tors in Eq. JSJ map onto three Hcrmitian 8x8 matri- 
ces obeying the same angular momentum commutators 
is given by the choice 



Vit) = (Pll-P33,-y|(/Oll+P33-2p22),Pl2+P21,P21 

-Pl2 , Pl3 + P31 , P31 - Pl3 , P23 + P32 , P32 - P23 ) • (4) 



Our choice differs only slightly from that in where 
this set is called a "coherence" vector. The resulting 
equation for ri{t) takes the form 



Z7)(<) = C{t)f^{t), 



(5) 



with 



£{t) = -in + e{t)B^ + 2J{t)B^, 



(6) 



where F indexes the strength of the decoherence. The 
matrices B take the form 



B, = 



/o 
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(7) 



As per the unitary integration procedure [3, [3j the 
solution of Eq. (|SJ is written as a product of exponentials 

r]{t) = exp[— Fi] exp[— 

X cxp[-i^_(t)S_] ex.p[-i^j,{t)B^]r){0), (8) 

with B± = Bx ± iBy. Because our procedure depends 
only on the commutation relations which remain as in 
2] , the classical functions /i in the exponents satisfy the 
same equations as before. 



ii+-ie{t)n+- J{t){\+ n\) 
/i = 2^J{t)^l+ - eit), 
fi- — ijiji^ = J(t), 



= 0, (9a) 
(9b) 

fi,{0) = 0. (9c) 



The first of these equations, involving p+{t) alone in Ric- 
cati form, is the only non-trivial member of this set. Once 
solved, /x_ and fj, are obtained through simple integration 
of the remaining two equations. For given e{t) and J{t), 
a mathematica program solves the set of equations read- 
ily. Also, the subsequent algebra involved in evaluating 
the exponentials in Eq. and their product is easily 
carried out. Thereby, for any initial density matrix and 
its ?7(0), we obtain rjit) and thus p(t) at any later time. 

Since our model for decoherence introduces its effect 
through the single real factor which is the first term on 
the right-hand side of Eq. ((HJ, the density matrix remains 
Hermitian throughout. This is an advantage, permitting 
evaluation of its eigenvalues and calculation of quanti- 
ties such as the entropy of the system. It is also clear 
that for any finite F all elements in ri{t) in Eq. vanish 
asymptotically with t so that all off-diagonal elements of 
the density matrix so vanish while all diagonal elements 
become equal. With the trace invariant and chosen to 
be unity, the density matrix evolves to that of the so- 
called chaotically mixed state, ^I. Correspondingly, the 
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entropy reaches asymptotically the value InS = 1.0986. 
These are aspects of the general result valid for all n- 
level systems . We note again that other models of de- 
coherence and dissipation through other choices for the 
operators Lfc in Eq. ^ than the one we made will, in gen- 
eral, lead to a larger set of exponential factors in Eq. ||SJ|, 
making for more complicated algebra therein and in the 
coupled set of equations in Eq. (9). However, inclusion 
of a term involving also Ay in Eq. (O, that is, a cou- 
pling also between levels 1 and 3, causes no additional 
difficulty since it does not enlarge the number of ^ or _B 
matrices in our procedure. 



TWO DIFFERENT DRIVING FIELDS BETWEEN 
NEIGHBORING STATES 

We present results for three degenerate states, such as 
of three identical pendula, with different nearest neighbor 
couplings between 1-2 and 2-3, that is, with e(t) and J(t) 
differing in amplitude and frequency, 



e{t) = Acosnt, J{t) ^ -Bcos{Lot + d), (10) 

with 6 a relative phase difference. A representative sam- 
ple of the density matrix upon starting with all popula- 
tion in the state 1 and all other elements zero is shown 
in Figs. 1-4. Note the appearance of a complicated fre- 
quency spectrum beyond just the two introduced driving 
frequencies. The analytically solvable problem presented 
in the next section provides an understanding of the ori- 
gin of these other frequencies. As shown in Fig. 3, the 
entire population can be transferred from level 1 to level 
3 over certain time intervals. 

Specializing to equal driving frequencies with a fixed 
amplitude ratio, results for various phase differences be- 
tween the two fields are shown in Figs. 5-8. Clearly, the 
density matrix elements depend on the relative phase. 

To contrast with a different initial state. Figs. 9-10 
show results when all population is in the state 2 at i = 0. 
Fig. 11 presents the evolution of the entropy, S — 
— Trplnp, showing a monotonic rise independent of am- 
plitudes and phases and of the initial pure state. In- 
deed, the eigenvalues of the density matrix are |(1 — 
e~'"*), ^(1 — e~'"*), and ^(l -I- 2e~'"*), from which the en- 
tropy easily follows. 



3 STATES OF THE HYDROGEN ATOM IN 
AN OSCILLATING ELECTRIC FIELD 



An example of a three-state degenerate system is pro- 
vided by the n = 3, m = states of the hydrogen atom. 
An oscillating electric field such as that of incident ra- 
diation couples 3s — 3p and 3p — 3d states, the dipole 
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FIG. 1: Time evolution of the elements of the density matrix 
of a n = 3 system driven by the fields in Eq. IjlOfl and Hamil- 
tonian in Eq. with fl = 0,lj ^ 1,5 = 0,A = 0.05, B = 
0.5, r = 0.02. Right hand panels show the off-diagonal ele- 
ments, two of which are imaginary and one real. 
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FIG. 2: same as in Fig. 1, except that A — 0.15. 



matrix elements being in the ratio a/2 : 1. Our results in 
this paper apply to this situation with the two frequen- 
cies m Eq. (Uni equal and A/B = V2. These were the 
choices made in Figs. 5-8. We present in Figs. 12-15 a 
sample of results for initial population in 3s for different 
amplitudes of the driving field. 

This problem is, of course, exactly solvable in terms 
of the parabolic eigenstates of hydrogen. With H{t) in 
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FIG. 3: same as in Fig. 1, except that A = 0.5, B = 1,0.= 
0.1,0; = 1. 
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FIG. 4: same as in Fig. 3 except that A= 1,B = l/\/2, F = 
0.08 and longer times shown to illustrate asymptotic evolu- 
tion. 



FIG. 5: Elements of the density matrix with driving fields 
of same frequency, fl = cj = 1, and amplitudes A = 1, B = 
I/V2, and r = 0.02, 5 = -7r/6. 
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Eq. I2I1 containing a single time dependence, the resulting 
Schrodinger equation, 



FIG. 6: same as in Fig. 5 except that 5 = tt/G. 



( m \ / -A \ / s{t) \ 

i pit) = -A/V2 p{t) coacjt, 

V dit) J V -A/V2 ; V dit) J 

(11) 

can be solved after diagonalizing the matrix of constant 
coefficients to obtain the parabolic eigenstates ± eigenstate is then easily followed 



j=d, -j={s — \/2d)} and corresponding eigenvalues 



A{± 4/1,0}. The independent time evolution of each 
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FIG. 7: same as in Fig. 5 except that 5 = tt/A. 
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FIG. 8: same as in Fig. 5 except that 5 = 7r/2. 



FIG. 9: Evolution of density matrix elements, starting with an 
initial non-zero value only for p22 ~ 1. Contrast with Figs. 
1-4. The parameters are Q = uj ^ 1, A ^ 2, B ^ V^,T = 
0.02,(5 = 0. 
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Thus, for initial population in the s state, we have 



1 /3 
s{t) = + 2cos[W-(A/w)sinwt]) (12a) 
o V ^ 

= ^ism[^{A/u;) sin u;t] (12b) 

(i(i) = ^(cos[A/^(A/cj)sincjt] - 1). (12c) 
o V ^ 

Together with the exponential decrease of the elements of 
7y(t), the density matrix can be constructed to reproduce 



FIG. 10: same as in Fig. 9 except that yl = 10, B = 5%/2. Note 
the appearance of more rapid oscillations with the harder 
driving fields. 



the results in Figs. 12-15. It is also clear that when one 
of the parabolic states is used as a starting point, the 
density matrix will remain frozen for F = and decay 
monotonically for finite F as shown in Figs. 16 and 17. 

A recent paper has presented results similar to the 
above in Eqs. (12) for n = 2,3 0. The results can be 
readily extended to any n since the expansion of parabolic 
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FIG. 11: Evolution of the entropy to accompany the results 
shown in previous figures. The rise is monotonic from to 
In 3, the rate of rise depending only on the value of F. 
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FIG. 13: same as in Fig. 12 except that A = 2. 
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FIG. 14: same as in Fig. 12 except that A = 10. 



FIG. 12: Evolution of the density of states of the n = 3 Stark 
field with initial population in the 3s state. The amplitude of 
the driving field is A = 1 and F = 0.02. 



states in terms of spherical states of hydrogen are well 
known and given by 3j-coefficients 10]. The occurrence 
of the "Floquet form" in Eqs. (12), with trigonometric 
functions whose arguments are themselves a trigonomet- 
ric function scaled by A/lu, accounts for the appearance 
of higher frequencies than uj in Figs. 1-8 for stronger driv- 
ing fields. 

We note that such studies of Rydberg atoms in an n 



manifold under microwave ionization, sometimes with an 
additional static field, have been of considerable experi- 
mental and theoretical interest llll Il3| . 



SUMMARY 

The method of unitary integration as extended to prob- 
lems involving dissipation and decoherence affords a con- 
venient and powerful way of treating n-state systems in 
time-dependent fields. Through the solution of a sin- 
gle, classical, Riccati equation (first-order in time and 
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FIG. 15: same as in Fig. 13 except with a larger V = 0.08. 
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FIG. 17: same as in Fig. 16 except that V = 0.2. Note the 
monotonic evolution from the initial Stark state to the mixed 
state described by the density matrix il. 
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FIG. 16: Evolution of the n = 3 states of hydrogen, starting 
with a Stark eigenstate, and F = 0. Such eigenstates remain 
frozen in their time dependence, in contrast to other initial 
states as shown in the previous figures. 



quadratically nonlinear), we can follow the evolution of 
the density matrix in time without any restrictions to 
infinitesimal stops or time ordcrings. Results have been 
presented for three-state systems with examples of cou- 
pled pendula and the n = 3 states of hydrogen in a radi- 
ation field. 
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